Enhanced dispersion interaction in confined geometry 
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The dispersion interaction between two point-like particles confined in a dielectric slab between 
two plates of another dielectric medium is studied within a continuum (Lifshitz) theory. The retarded 
(Casimir-Polder) interaction at large interparticle distances is found to be strongly enhanced as the 
mismatch between the dielectric permittivities of the two media is increased. The large-distance 
interaction is multiplied due to confinement by a factor of (337 5//2 + 137 _3//2 )/46 at zero temperature, 
and by (57 2 + 7~ 2 )/6 at finite temperature, 7 = ei n (0)/e ou t(0) being the ratio between the static 
dielectric permittivities of the inner and outer media. This confinement-induced amplification of 
the dispersion interaction can reach several orders of magnitude. 
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The dispersion interaction acts between any two po- 
larizable objects, thus being one of the most ubiquitous 
interactions in nature 0, 0. It plays a central role in 
numerous phenomena in chemical physics and materials 
science, including gas-liquid condensation, capillarity 
inter-surface interactions [4J, and colloid stability [5j- 

The dispersion interaction is a quantum fluctuation- 
induced coupling between two polarizable particles me- 
diated by the electromagnetic field. In an unconfined 
system at zero temperature there is a single length scale 
with which the interparticle distance R is to be com- 
pared, i.e., the characteristic wavelength Ao of photon 
absorption by the particles, typically in the ultravio- 
let to visible range. London's calculation 0, valid 
in the nonretarded limit R <C Ao, yields the poten- 
tial U(R) = -[(3K/ir) r o °° d£a 2 {i£)]R- 6 , where a(u) is 
the frequency-dependent polarizability of the particles. 
Casimir and Polder recast the problem in quantum- 
electrodynamic terms, whereby the interaction arises 
from the effect of the particles on the zero-point modes 
of the electromagnetic field. For R <C Ao the Casimir- 
Polder result coincides with London's, yet in the re- 
tarded limit, R 3> Ao, the interaction decays as i?~ 7 , 
U(R) = -[(23/(47r))?ica 2 (0)]ir 7 8]. As a result, the 
dispersion interaction between two particles in the re- 
tarded regime (typically R > 0.1 /xm) is extremely weak 
and has not been directly observed. (Particle-surface and 
surface-surface interactions across such micron-scale dis- 
tances are much stronger and were successfully measured 
in the 1990s 0, Hc|.) In the current Letter we demon- 
strate that this weak particle-particle interaction can be 
dramatically amplified in confined geometries. 

At a finite temperature T another length scale appears, 
i.e., the thermal wavelength Xt — Tic/T ~ 7.6 /xm at 
room temperature. (The Boltzmann constant is set here- 
after to unity.) As was shown in Ref. ^H, retardation and 
finite-temperature effects arc intertwined. For R At, 
the Helmholtz free energy of interaction returns to a R~ 6 



dependence, F(R) = -3Ta 2 (0)R- 6 . 

The theory of dispersion interactions was extended by 
Lifshitz et al. to the case where the interacting objects 
as well as the intervening space are continuous media 
|l2l Il3| . In this continuum theory the material response 
to electromagnetic fields is assumed to be fully captured 
by the com plex , frequency-dependent dielectric permit- 
tivity e(<x>) |l4j . The current work is based on such a 
continuum approach. Despite the strong underlying as- 
sumption (after all, the media themselves consist of dis- 
crete polarizable particles), the Lifshitz theory has been 
widely used and experimentally corroborated 0. It is 
expected to yield valid results so long as the distance 
between the two particles is much larger than the inter- 
molecular distances in the materials. 

In various circumstances particles are spatially con- 
fined, e.g., in porous media, micro-cavities, biological 
constrictions or nanofluidic devices. Such confinement in- 
troduces a new length scale, the separation h between the 
bounding surfaces. The dispersion interaction between 
a single particle and confining surfaces has been exten- 
sively studied in the context of cavity QED 0, E>] • The 
effect of confinement on the interaction between two par- 
ticles, however, has been only partially addressed. Con- 
finement by two metallic plates (i.e., the limit where the 
permittivity of the outer medium e out — * 00) was found to 
drastically affect the interaction between two point-like 
particles in vacuum at zero |l6T ] and nonzero |l7| temper- 
atures. The nonretarded interaction between particles 
confined by two dielectric plates at T = was addressed 
within a single- image approximation in Ref. |18| . We fol- 
low the lines of Mahanty, Ninham, Bostrom and Longdell 
[lil Kl\ and extend their theory to the general and more 
practical case of arbitrary permittivities of both inner 
and outer media, discovering a dramatic enhancement 
for experimentally relevant values of e ln and e ou t. 

The system under consideration is schematically shown 
in Fig. ^ Two point-like, isotropic particles of polariz- 
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ability a(uo) are embedded in a slab of thickness h and 
dielectric permittivity ei n (w). The slab is bounded by two 
semi-infinite media of dielectric permittivity e ou t(o->). We 
use cylindrical coordinates, r = (p, ip, z), the z axis taken 
perpendicular to the bounding surfaces. For simplicity 
we specialize to the symmetric case where the particles lie 
on the slab midplane, z — h/2, connected by the vector 
R = (R, 0,0). All materials are assumed nonmagnetic. 



r = (R, 0, h/2), the tensor becomes diagonal, 



"out 



"m 



R 



z=h 



z=0 



"out 



FIG. 1: Schematic view of the system and its parameters. 

We employ the semiclassical scheme introduced by Ma- 
hanty and Ninham, which accurately reproduces the Lon- 
don and Casimir-Polder results 19] . In this theory the 
energy of interaction at T = is given by 

i-oo m 

T = : U(R) = -8-Kh / d^a 2 (i£) 
xTr[G(R,iOG(-R,ie)], (1) 

and the Helmholtz free energy at finite T by 

oo 

T > : F(R) = -(^T^VfeX^/c 4 ) 

xTr[G(R,i£„)G(-R,i£„)L = (2*c/\ T )n, (2) 

the prime indicating that the n = term is multiplied by 
1/2. In Eqs. JTJ and © R = r - r', and G(r,r',w) is the 
dyadic Green tensor of the electric-field wave equation, 



V x V x G - e m (oj)(uj 2 /c 2 )G = IS(r - r'), 



(3) 



where I is the identity tensor and e m = €i n or e out de- 
pending on whether r lies in the inner or outer medium. 
(The position r' is taken inside the slab.) The boundary 
conditions in the current case are continuity across the 
bounding surfaces of the tangential components of both 
the electric and magnetic fields. This imposes continuity 
on z x G and z x V x G across z = and z = h. 

Thus, this scheme reduces the problem to finding the 
Green tensor G(r,r',uj) of Eq. © with the aforemen- 
tioned boundary conditions. The derivation is techni- 
cally complicated and can be found in Refs. [U |H|. 
In the symmetric case of interest, r' = (0,0, h/2) and 



dj (R, i£) = / dxg l (R, if, x)<% 

„2u 



g p = [q 2 x(e v - t){e v - u)J (x) - (x\e lv + tu) 
-eV + q 2 )(t + u))Mx)]/[q{e" - t)(e v + «)] 

g v = [p 2 x{e v + t)(e v + u)J a (x) + (x 2 {e 2v + tu) 
-e V + q 2 )(t + u)) Mx)]/[q(e v - t)(e- + «)] 

g z = -x 3 (e v +u)J (x)/[q(e v -u)}, (4) 

where Jk(x) are Bessel functions, and the following ab- 

1/2 

breviations have been used: 7 = EmAout, P = R£ e in l c > 
q = {x 2 +p 2 )^ 2 , s = (x 2 + p 2 /i) l/2 , t= (q-s)/(q+s), 
u = (q - js)/(q + 7 s), and v = hq/R. Given e in , ou t(if) 
and a(if), one can substitute Eq. J3J in Eq. Q or (0) 
and calculate numerically the interaction potential. 

It is instructive, however, to first analyze the inter- 
action in several asymptotic limits. We begin with the 
small-distance limit, recovering the known results for 
unconfined particles. For R <C h the expressions for 
gi are expanded to leading order in large v, where- 
upon the integration in Eq. (@J can be carried out an- 
alytically. Substituting the result in Eq. QJ, we get 
U(R) = ~(h/ir)R- e r °°^e- 2 P(3 + 6p + bp 2 + 2p 3 + 
p 4 )a 2 (if)/e 2 n (if). In the nonretarded limit, R -C A , we 
take the leading order in small p, for which the London 
result is recovered, 



T = 0, R < h, A : U(R) = 



3h [°° a 2 ^) 

with the appropriate correction due to the fact that the 
particles are not in vacuum (e ul 7^ 1). In the retarded 
limit, R 3> Ao, the integral over £ is dominated by small 
(static) £. An asymptotic analysis then recovers the 
Casimir-Polder result with the correction for ei n 7^ 1, 



T = 0, A < R < h : U(R) 



23^ca 2 (0) 
^ef n /2 (0)i? 7 



(6) 



Similarly, substitution of the asymptotic 
Green tensor for R <C h in Eq. © yields 
F(R) = -2TR- 6 J2'™ e- 2 ^(3 + 6 Pn + 5p 2 + 2p\ + 

pi)a 2 (iCn)/e 2 n (i£n), where p n = R£ n e\l 2 /c. Subse- 
quently taking the limit R <C At, we expand to leading 
order in small p n |22l | and recover the known result for 
the nonretarded interaction at finite temperature |2.'"{l| . 



T > 0, R < h, A T : F(R) 



6T 



a 2 (if„ 



^0 



(7) 



In the other limit, i? 3> At, the sum is dominated by the 
n = term, leadin g to the known retarded interaction at 
finite temperature |ll|. 



T > 0, A T < R < /i : 



3ra 2 (0) 
£ 2 n (0)i?6- 



(8) 
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We now turn to the more interesting asymptotic limit 
of large R where confinement sets in. In the limit R^> h 
the expressions for gi in Eq. 10} are expanded to lead- 
ing order in small v, and the integration in Eq. Q 
is performed analytically. The result is substituted in 
Eq. □ to yield U(R) = -(h/(2ir))R- 6 d£ e - 2 P[(5 + 
7~ 4 ) + (10 + 2 7 - 4 )p + (7 + 3 7 - 4 )p 2 + (2 + 2 7 ~ 4 )p 3 + (1 + 
7~ 4 )p 4 ] a2 (*£)/ e out(*£)j where p = 7 ~ 1 / 2 p. In the non- 
retarded regime, this expression is expanded to leading 
order in small p, leading to 



T = 0, /i< R < A : 

3h 



(9) 



U(R) = - 



irR 6 



5 7 2 ( ? o+ 7 - 2 ^) 



Equation JJj} gives the nonretarded interaction under 
confinement. It reduces to the London result, Eq. J5J, 
for 7 = 6i n /e out = 1. This result has a rather restricted 
validity as it requires that h be much smaller than Ao- 
(A detailed discussion of the confined nonretarded inter- 
action will be given elsewhere.) In the retarded regime 
the integral is dominated by small £, leading to 



T = 0, R > h, A : 

23W(0) 33 7 5/2 



13 7 



-3/2 



47ref n /2 (0)^ 



40 



(10) 



where we have written 7 (0) = 7 for brevity. This re- 
tarded interaction under confinement converges to the 
Casimir-Polder expression, Eq. ©, for 7=1. 

At finite temperature we substitute the asymptotic 
Green tensor for R ^> h in Eq. (0), resulting in F(R) = 

-TR- 6 EZo^ 2pn i^ + 7~ 4 ) + (10 + 2 7 - 4 )p„ + (7 + 

3 7 - 4 )p 2 + (2 + 2 7 - 4 )^ + (1 + l~ 4 )P 4 n ]a 2 (ii; n )/t 2 ouMn), 
1/2 

where p n — R^ n eJ^lc. For R <C At we take the leading 
order in small p n [22j and get 

T>0, /i<i?< A T : (11) 
6T W a 2 (z£„) 5 7 2 ^„) + 7~ 2 fe) 



n=0 
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This is the extension of Eq. to the confined case. For 
R 3> At, the n = term dominates the sum, yielding 



T > 0, i? > h,X T : F(R) = 



3Ta 2 (0) 5 7 2 + 7- 



MR 6 



6 



(12) 

which extends Eq. JHJ to the confined geometry. 

Equations (fTUfl and ifT^ are our central results. They 
account for the large-distance, retarded dispersion inter- 
action between the confined particles at zero and finite 
temperature, respectively. Comparing with Eqs. © and 
(JBJ, we see that the confinement is manifest as a factor 
dependent on the ratio 7(0) between the static dielectric 
permittivities of the inner and outer media. This factor 



can be as small as 0.78 [Eq. (JTOJl] or 0.75 [Eq. |0] , but in- 
creases indefinitely with the permittivity mismatch. The 
divergence of the interaction energy for 7 — » or 00 
is obviously unphysical. Although Eqs. I|10l) and l|12|) 
are asymptotically correct for any finite mismatch, as 
7 becomes increasingly large or small one must go to 
ever larger interparticle distances for these asymptotes 
to hold. Ultimately, in the limits 7 — > 0, 00 their range 
of validity disappears, and the large- distance interaction 
obeys a different power law 

[HEi (Detailed analy- 
sis of this behavior will be given elsewhere.) The main 
point, however, is that the amplification factor can be 
very large for reasonable values of 7. For example, for 
particles embedded in a polar liquid [ei n (0) = 80] which 
is confined by glass plates [e ou t(0) = 4] at room temper- 
ature, one gets an amplification factor of about 300. If 
the outer medium is a gas (a free-standing film, e out = 1), 
the factor increases to about 5000 [24| . 

Finally, we present results from numerical integration 
of Eqs. @ and Q for two examples of practical inter- 
est. This requires expressions for e(i£) of the various 
media, for which we use the Ninham-Parsegian repre- 
sentation 0, [2f| , an empirical fit based on electromag- 
netic absorption spectra of the materials. In the first 
example two polystyrene particles are confined in a slab 
of water between two glass plates at room temperature. 
The function ei n (^) for water is found in Refs. 0,0 and 
that for silica glass, e out (i£,), in Ref. [2j|. For a(i£) we 
took the excess Clausius-Mossotti polarizability |j] of a 
polystyrene sphere, a(iQ/V = [3e in (^)/(47r)][eps(^) - 
e in(*£)]/[ e Ps(j£) + 2ei n (i£)], V being the particle volume 
and eps(*£) the permittivity of polystyrene, found in 
Refs. |2j, [5| . The resulting potentials for two interplate 
separations, h = 1 and 0.1 /im, are shown in Fig. EJa) 
along with the unconfined potential. The interaction per 
volume squared has been scaled by —T/R e . The curves, 
therefore, represent the effective Hamaker coefficient (di- 
vided by 7r 2 , in units of T) as a function of interparticle 
distance. The unconfined potential clearly exhibits the 
crossover from the nonretarded i?~ 6 regime at small R to 
the retarded R~ 7 dependence at intermediate distances, 
and then back to the R~ 6 decay due to temperature. 
The confined interaction deviates from the unconfined 
one at distances R ^ h. The amplification factor in- 
creases moderately with distance until saturating to the 
asymptotic value given in Eq. (|12fl . In the second ex- 
ample two polystyrene particles are confined in a hydro- 
carbon slab of thickness h = 3 nm, which is embedded 
in water. This may mimic small hydrophobic inclusions 
in a biological membrane. For the oily environment we 
used the permittivity £i n (*£) of pentane Q. As is seen 
in Fig. [3(b) , the amplification becomes significant only 
at R > lOh. This is because for smaller distances retar- 
dation has not yet set in. (Note where the unconfined 
potential departs from its nonretarded R~ 6 behavior.) 

The confinement-induced enhancement of the retarded 
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FIG. 2: Potentials of interaction as obtained from numerical 
integration of Eqs. (2) and (|1J. The interaction free energy 
is scaled by —T/R 6 . (a) Polystyrene particles in water be- 
tween two glass plates at T = 300K. Solid, dashed, and dot- 
ted curves correspond, respectively, to an unconfined system, 
h — 1 fim, and ft = 0.1 fj.m. (b) Polystyrene particles in a 
hydrocarbon slab embedded in water at T = 300K. Solid and 
dashed curves correspond to an unconfined system and h = 3 
nm, respectively. Insets in both panels present the ratio be- 
tween the confined and unconfined potentials, the dash-dotted 
lines showing the asymptotic amplification factor of Eq. 1121 1. 

dispersion interaction can be viewed as a consequence 
of multiple-reflection waveguiding of the electromagnetic 
radiation between the boundaries. Despite the demon- 
strated strong effect the interaction remains weak. In 
the example presented in Fig. |2Ja), for instance, the in- 
teraction free energy for h = 1 /im and R = 10 /xm 
is about 3 x 10~ 6 T//j,m 6 . Nevertheless, to achieve the 
same energy without confinement one would have to set 
the interparticle distance at about 2 /zm. Thus, besides 
the fundamental significance of the strong confinement 
effect reported here for particle interactions in confined 
systems, it may become useful also in extending the range 
of observation of the Casimir-Polder interaction. 
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